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Interaction of the internal states of an atom with spatially dependent standing-wave cavity field 
can impart position information of the atom passing through it leading to subwavelength atom lo- 
calization. We recently demonstrated a new regime of atom localization [Sahrai et ai, Phys. Rev. 
A 72, 013820 (2005)], namely sub-half- wavelength localization through phase control of electromag- 
netically induced transparency. This regime corresponds to extreme localization of atoms within a 
chosen half-wavelength region of the standing-wave cavity field. Here we present further investiga- 
tion of the simplified model considered earlier and show interesting features of the proposal. We show 
how the model can be used to simulate variety of energy level schemes. Furthermore, the dressed- 
state analysis is employed to explain the emergence and suppression of the localization peaks, and 
the peak positions and widths. The range of parameters for obtaining clean sub-half-wavelength 
localization is identified. 

PACS numbers: 42.50.Ct, 42.50.Pq, 42.50.Gy, 32.80.Lg 
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I. INTRODUCTION 

Precise localization of atoms has attracted consider- 
able attention in recent years. Optical manipulations 
allow probing the center-of-mass degrees of freedom of 
atoms with subwavelength precision. The interest in sub- 
wavelength atom localization is largely due to its appli- 
cations to many areas requiring manipulations of atomic 
center-of-mass degrees of freedom, such as laser cool- 
ing jBose- Einstein condensation and atom lithog- 
raphy 3] alongwith fundamentally important issues such 
as measurement of the center-of-mass wavefunction of 
moving atoms 

Optical techniques for position measurements of the 
atom are of considerable interest from both theoretical 
and experimental point of view. Several scheme have 
been proposed for the localization of an atom using op- 
tical methods p|. It is well known that optical methods 
provide better spatial resolution in position measurement 
of the atom. For example, in the optical virtual slits 
scheme the atom interacts with a standing-wave field 
and imparts a phase shift to the field. Measurement of 
this phase shift then gives the position information of the 
atom Q . In another related idea based on phase quadra- 
ture measurement is considered in Ref. |7(. Kunze et 
al. @ demonstrated how the entanglement between the 
atomic position and its internal state allows one to local- 
ize the atom without directly affecting its spatial wave 
function. It is shown that, by using Ramsey interfer- 
ometry, the use of a coherent-state cavity field is better 
than the classical field to get a higher resolution in po- 
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sition information of the atom Resonance imaging 
methods have also been employed in experimental stud- 
ies of the precision position measurement of the moving 
atoms 00. 

More recently atom-localization methods based on the 
detection of the spontaneously emitted photon during 
the interaction of an atom with the classical standing- 
wave field, are considered 0,0,0,0- It is, however, 
important to note that from an experimental point of 
view, observation of spontaneous emission spectrum is 
very tricky and difficult. In this context, another scheme 
based on a three-level A-type system interacting with two 
fields, a probe laser field and a classical standing wave 
coupling field, is used for atom localization by Paspalakis 
and Knight [lfj . They observe that in the case of a weak 
probe field, measurement of the population in the up- 
per level leads to sub- wavelength localization of the atom 
during its motion in the standing wave. Thus, in essence, 
this scheme uses absorption of a probe field for atom lo- 
calization. Atomic coherence effects, such as coherent 
population trapping, have also been shown to be useful 
for subwavelength localization of atoms by Agarwal and 
Kapale ^3] i where monitoring the coherence of the trap- 
ping state gives rise to subwavelength localization of an 
atom to a precision prechosen through the ratio of the 
square of Rabi frequencies of the strong standing-wave 
drive field and a weak probe field. 

The authors (with collaborators) recently proposed a 
subwavelength atom localization scheme through phase 
control of the absorption of a weak probe field by the 
atom. A modified A-type level scheme with an extra 
level and the drive fields forming a complete loop was 
shown to introduce a phase dependence in the response 
of the atomic medium to a weak probe field. This phase 
controllable atomic response was shown to give rise to 
tunable group velocity from subluminal to superluminal 
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in a single system 18] . By considering one of the drive 
fields to be a standing- wave field of a cavity it was shown 
that the same scheme can be used to localize an atom 
flying through the standing wave field to subwavelength 
domain [19|. 

This article is a sequel to the earlier article jl9( , hence- 
forth referred to as J* . In J> a restricted parameter range 
of the model was considered to show the possibility of 
sub-half-wavelength localization. In this article further 
investigations of the analytical results are carried out to 
show how the model can be used to simulate variety of 
atomic systems with varying energy level spacings, dif- 
ferent atomic dipole matrix elements and decay proper- 
ties. Appropriate parameters required to obtain different 
regimes of localization are studied in detail. A dressed- 
states approach is also considered to give insight into the 
results obtained. 

The article is organized as follows: For completeness, 
a brief description of the procedure to determine the sus- 
ceptibility of the atom to a weak probe field is given in 
Sec.[n] Then, in Sec. IIIII A. the susceptibility expression 
is studied in detail to arrive at the conditions for ob- 
serving atom localization. Various parameter ranges for 
the drive field Rabi frequencies and decay properties are 
considered in order to simulate variety of atomic species 
and to clarify experimentally controllable features and 
properties of the model and numerical results and their 
explanation through the analytical probe susceptibility 
expression is presented in Sec. 11111 B. A simple dressed 
states treatment is presented in Sec. IIIIC in order to ex- 
plain the results obtained in the earlier sections. Finally 
the conclusion is presented. 
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FIG. 1: The Model: (a) The cavity supports the standing 
wave field (1) corresponding to Rabi frequency Oi. Two other 
fields (2, 3) are applied at an angle as shown. The atom enters 
the cavity along the z axis and interacts with the three drive 
fields. The whole process takes place in the x-z plane, (b) The 
energy level structure of the atom. Probe field, denoted by 
£ p , is detuned by an amount A from the |ai) — \c) transition. 
The fields (2, 3) shown in (a) part of the figure correspond to 
the fields with Rabi frequencies fi2 and £^3 respectively. The 
decay rates from the upper levels |ai) and I02) are taken to 
be 71 and 72 respectively. 



bution of the atom is nearly uniform along the direction 
of the standing wave. Therefore, we apply the Raman- 
Nath approximation and neglect the kinetic part of the 
atom from the Hamiltonian 20]. Under these circum- 
stances, the Hamiltonian of the system in the rotating 
wave approximation can be written as 



(1) 



where 



Ho = tuJ ai \ai) (ai\+huj a2 \a 2 ) (a 2 \+Tuo b \b) (b\+huj c |c) (c| 

(2) 

and 



II. THE MODEL AND EQUATIONS 

The schematics of the proposed scheme are shown in 
Fig. ^ We consider an atom, moving in the z direction, 
as it passes through a classical standing-wave field of a 
cavity. The cavity is taken to be aligned along the x axis. 
The internal energy level structure of the atom is shown 
in Fig.^b). The radiative decay rates from the level |ai) 
and 1 a 2 ) to level |c) are taken to be 71 and 72. The upper 
level \a\) is coupled to the level \a%) and further the level 
1 02) is coupled to level \b) via classical fields with Rabi fre- 
quencies O3 and r?2, respectively. In addition, the upper 
level |ai) is coupled to level \b) via a classical standing- 
wave field having Rabi frequency fii. It should be noted 
that the Rabi frequency of the standing wave is position 
dependent and is taken to be Cli(x) = Jlisinra . Here 
Qi(x) is defined to include the position dependence and 
k is the wave vector of the standing wave field, defined 
as Hi = 2tt/X, where A is the wavelength of the standing- 
wave field of the cavity. We assume that the atom is 
initially in the state |c) and interacts with a weak probe 
field that is near resonant with |c) — > |oi) transition. The 
detuning of the probe field on this transition is taken to 
be A. We assume that the center-of-mass position distri- 
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Here u>i are the frequencies of the states |z) and i/j are the 
frequencies of the optical fields, and 9 2 , 63 are the angles 
made by the propagation direction of the fields £l 2 and 
Q3 with respect the x axis respectively. The subscript 
p stands for the quantities corresponding to the probe 
field — i.e., £ p and v p are the amplitude and frequency of 
the probe field. Also p QlC is the dipole matrix element 
of the \c) — > \a±) transition. For simplicity, we assume 
that the Rabi frequencies f^i and Sl 2 are real and is 
complex — i.e., = Q.3e~ lv . This choice of imparting a 
carrying phase to field 3, is only for the convenience of 
calculations. As will become clear later, only the relative 
phase of the three fields is important and absolute phases 
do not matter. The dynamics of the system is described 
using density matrix approach as: 
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-[H,p]--{T,p} 



(4) 
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where {T, p] = Tp + pT. Here the decay rate is incorpo- 
rated into the equation by a relaxation matrix T, which is 
defined, by the equation (n|r|m) = 7„<5„ m . The detailed 
calculations of these equations are given in the Appendix 
oiS. 

Our goal is to obtain information about the atomic po- 
sition from the susceptibility of the system at the probe 
frequency. Therefore, we need to determine the steady 
state value of the off-diagonal the density matrix element, 
p ai c- After necessary algebraic calculation and moving to 
appropriate rotating frames, we obtain a set of density 
matrix equations. To determine p aiC we only need fol- 
lowing equations 

Pa lC = -[i(w„ ic - v p ) + \ii]p ai c + ^ 3 e-^e ikxcosf>3 p a2C 

+ -ill SUXKXpbc ~ 1 2 ^ (Pa iai - Pec), 
Pa 2 c = -[i(Wo a c ~ (Vp - V3)) + lp2\Pa 2 c 

-1— Poaai, 

Pbc = -[i(w&c + 1/1- i/ p ) + 7hc]pbc + sin«;a;p aic 



+ 2 SZ 2e Pa 2 c-1 2 ^ P6ai- 



(5) 



As we know, the dispersion and absorption are related 
to the susceptibility of the system and is determined by 
Paxc- We take the probe field to be weak, and calculate 
the polarization of the system to lowest order in £ p . We 
keep all the terms of the driving fields but keep only 
linear terms in the probe field. The atom is initially in 
the ground state |c), therefore we use 

5(0) _ 1 =(o) _ n ~(o) _ n ~(o) _ n (c] 

Equation (JSj can then be simplified considerably to ob- 
tain 
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itp ikx cos #3 ~ 



Pa 2 



Pa lC = -(iA + -ji)p aiC + e lip e 

i ~ . EpPaxc 

+ -ili sin Kxp bc + 1 2 ^ , 
/5 a2C - -(iA + i 7 2)/5 a2C + ^ 3 e iv e- itecose ^ QlC 



p bc = -i A pbe + -fii sin p QlC 



„"2 e Pa 2 c- 



(7) 



Here we have introduced the detuning of the probe field 
and the frequency difference between levels \ai) and |c), 

A = UJ aiC - V p = L0 a2C + Vz~V p = Ubc +v± — v p . (8) 



Here we have also assumed that -fbc = 0. It can be eas- 
ily seen that these set of equations can also be used to 
simulate a variety of level schemes as shown in Fig. El 
after redefining the decay rates accordingly as discussed 
in the caption. The schme in Fig. |2f&) requires special 
attention as the positions of the states |ai) and |a2) are 
reversed in the order of increasing energy compared to 
the other levelschemes. This entails small change in the 
rotating frame that is chosen to arrive at the simplified 
density matrix equations. The transformation required 
can be accomplished by replacing the complex Rabi fre- 
quency r?3 by its complex conjugate fil and changing its 
frequency ^3 to —1/3. The density matrix equations so 
obtained are identical to the set Q) given above except 
for the redifinition of the phase from tp — > —ip. How- 
ever, as will be seen later, the phase enters through the 
term cos <p in the response of the atoms to a weak probe 
field, thus the final results are identical for all the models 
discussed in Fig. [21 

This set of equations can be solved analytically; the 
detailed discussion can be found in the appendix of 
.y . Thus, the off-diagonal density-matrix element cor- 
responding to the probe transition is obtained as 



Pa lC = Pa lC e^ pt = —(nl- 4A 2 + 2i 72 A)£ p p aiC e- 

where we have chosen, without loss of generality, (73 
7r/4, 62 — 7r/2 + 7r/4, moreover, Y is defined to be 



(9) 



Y = A + iB, (10) 



with 



A = -8A 3 + 2A(ft 2 sin 2 kx + + nj) 

+ 2 7 i 72 A + nin 2 Q 3 (e iip + e~ iip ) sin kx, 
B = 4A 2 ( 7l + 72 ) - 7l f7 2 - 72 fi 2 sin 2 kx. (11) 

The susceptibility at the probe frequency can be writ- 
ten as 



_ 2Np aiC p aiC ^ t 



2N\p aiC \ 2 (fl 2 2 - 4A 2 + 2i 72 A) 



Yh 



(12) 

where N is the atom number density in the medium. The 
real and imaginary parts of susceptibility are given as 

x' = m e ^f m ~ AA2)A + 272AB} ' (13) 

" 2 ^%^{27 2 A^ - (n 2 - 4A 2 )B)}, (14) 



X 



where Z = YY* and \ — x' + l x" ■ It i s imperative to 
point out that the phase enters the susceptibility expres- 
sion only through the quantities A and Y . Even the phase 
dependence of Y is only through the quantity A. More- 
over, we observe that the phase dependent term in A is 
ri 1 £l 2 ^3(e lv3 + e -lv )sinra. Thus the phase factor could 
very well have come from either of the three driving fields. 



4 



(a) (6) (c) (d) 




FIG. 2: Several level schemes that can be studied using our model in Fig. 1(6), within the weak probe limit. The decay 
rates are defined as follows: (a) 71 = -y ai b + 7a iC and 72 = j a2 b + Ja 2 c; (b) 71 = j ai b + 7a iC and 72 = r y a2 a l + f a2 b + 7a 2C ; (c) 
71 = 7aia 2 + 7aib + Ja ic and 72 = 0; (d) 71 = 7 ai a 2 + 7aib + 7aic and 72 = 7a 2 &. Here 7^ corresponds to spontaneous decay 
rate from level \i) to level It can be noted that the positions of levels |ai) and |az) in (b) are reversed compared to the 
other schemes. Slight modifications in the equations are needed to simulate level scheme in (6) with the equations given in the 
text. The transformations required is: the complex Rabi frequency ^3 — > — i.e., vz — > —^3. It can be easily shown that the 
results remain unchanged under these transformations, as discussed in the text. 



As pointed out earlier, if all the fields had phase depen- 
dence, only the collective phase would be important and 
no individual phase-dependent terms would occur. This 
is because the Rabi frequencies fli in all the other terms 
appear through fi|, which is | Qi\ 2 for a complex Rabi 
frequency fli — |Qi|e l ^ >i . The collective phase can be eas- 
ily determined to be tp = <p2 + 03 — <f>i, by repeating 
the susceptibility calculation. Here 4>i is the phase of the 
complex Rabi frequency Qj of the ith driving field. 

In the next section we consider the imaginary part of 
the susceptibility x" m detail and obtain various condi- 
tions for subwavelength localization of the atom. 

III. RESULTS AND DISCUSSIONS 

We study the expression l|14|) for the imaginary part of 
the susceptibility on the probe transition in greater de- 
tail in the following discussion. It is clear that x" — i- e -> 
probe absorption — depends on the controllable parame- 
ters of the system like probe field detuning, amplitudes 



and phases of the driving fields. First we present analyt- 
ical considerations of the probe absorption maxima and 
its relation to the atom localization. Then we present 
the results of the numerical study for a variety of dif- 
ferent sets of values of the parameters. In the end we 
present the dressed-state analysis to shed some light on 
the numerical results. 



A. The probe absorption maxima 

Noting the dependence of x" on sinra, it is, in princi- 
ple, possible to obtain information about the x position 
of the atom as it passes through the cavity by measuring 
the probe absorption. Nevertheless, for precise localiza- 
tion of the atom the susceptibility should show maxima 
or peaks at certain x positions. We obtain the condi- 
tions for the presence of peaks in x" m the discussion 
to follow. Eq. (|14f> can be rewritten as follows, using 
,yV = 2N\p aiC \ 2 /he , 



X " _ A + B(kx) 

3r ~ 7l [A + 2B{nx)\ + 7 | (Q 2 sin 2 nx - 4A 2 ) 2 + [8A 3 - 2 A (fl 2 sin 2 KxTofTo|) - 2f2ift 2 f2 3 costp sinra] 2 

A + B(kx) 



71 [A + 2B(kx)} + -y 2 flf(smnx - Ri) 2 (smnx - R 2 ) 2 + 4A 2 S!f (sin kx - i? 3 ) 2 (sin - i? 4 ) 2 

where 



(15) 



.4 


= 71 (4A 2 - 


B(kx) 


= 72(^1 &l 




2A0 3 / 


L\.2 




Rl,2 


2A 





B(kx) = 72(^1 ^2 sm2 KX + 4 A fli f2 2 ^3 cos tp sinra + 4 A 2 Q 2 ) = 72 fl 2 £l\ (sinra — Li)(shiKx — L 2 ) , 

(- cos tp ± \J cos 2 tp - lj , (16) 



^3,4 = 



1 

2A0! 



-n 2 si 3 cos ip ± a/s^i st. 2 cos 2 99 - 4A 2 (si 2 + si 2 - 4A 2 ) 



It can be seen that the probe field absorption would peak 
at positions satisfying 

Sin KX = i?l,2,3,4 ■ (17) 

The roots Li,2 do not contribute to the probe absorption 
maxima as they are the roots of the numerator as well, 
and these contributions mutually cancel. Moreover, for a 
given set of parameters not all four roots contribute to the 
probe absorption maxima, as they have different weight- 
ing factors given by 72^1 and 4A 2 Sli. The dominant 
weighting factor, being independent of position, governs 
which set of roots {-^1,2} or {i?3,4} will be important for 
dictating the atom- localization positions. 



It can be clearly seen that for 72 = the maxima po- 
sitions gverned by R12 do not occur; whereas, for A = 
the roots i?i j2 are more important compared to i?3,4. In 
the regime where both A and 72 are zero interesting con- 
sequences follow. This competition of the roots gives rise 
to various interesting regimes of parameters and possi- 
bilities in the atom localization. In the following we will 
throw light on the novel propreties arising due to this 
freedom. It can be noted that in J 2 " the roots L12 and 
i?i,2 did not appear as 72 was taken to be zero, which 
leads to B(kx) — 0. For completeness we give the ex- 
pression of x" as use d in J^: 



2N\pa lC \ 2 

he 7i(0| 

2N\pa lC \ 2 



4A 2 



7i ("I 

(8A 3 - 2 A (ST. 2 sin 2 KX ' 
7i(^-4A 2 ) 2 



4A 2 ) 



2\2 



SI 2 ) — 2 f7iSl2^3 cosy shiKrr) 2 



he Q 7 2 (S1 2 -4A 2 ) 2 + 4A 2 Slf (sinra - i? 3 ) 2 (sin/ 



(18) 



A direct calculation of x" from the equation (|14fl shows 
that the positions of maxima of do not strongly depend 
on the decay parameters, and are function of only the 
drive field Rabi frequencies and phases. However, for 
a chosen value of the detuning the widths of the peaks 
observed in the plots of x" vs KX depend on the values 
of the decay parameters. To make connection with the 
positions of maxima predicted by the roots of the denom- 
inator in Eq. (|15fl . i.e, the roots i?i,2,3.4, we study these 
roots in detail in the following discussion. 

The probe field detunings required to obtain probe 
field absorption peaks as a function of the x coordinate 
along the cavity field axis can be obtained by solving 
equations sinra = i?i,2,3,4 for A. We denote the so- 
lutions for s'mKx = i?i,2 as 5% 2 and the solutions for 
sinKx = i?3,4 as £3,4,5. It can be easily shown that 



0i,2 = T - ^- sin kx 



(19) 



and #3,4,5 are the solutions of smnx = i?3,4 — i.e., the 
cubic equation, 

4# 3 - #(fi 2 sin 2 kx + fi 2 , + SI 2 ) - S7.if2 2 ^3 sin kx cos ip = 0. 

(20) 

When the relative phase cp = 7r/2, the above equation 
can be readily solved to give 



# 3 = o, #4,5 =±£ \/njsm 2 Kx + ni + nj. (21) 



Thus, for ip = 7r/2 the above equations give the val- 
ues of the probe detuning for observing probe absorp- 
tion maxima as a function of the spatial position along 
the standing-wave field. It is clear that for 6 = A = 
there is no atom localization possible as the probe ab- 
sorption would be the same at all spatial positions. We 
do not give expressions for #3,4,5 for the case of ip = Oi as 
they are quite complicated; however, they can be readily 
evaluated numerically to verify the predictions. 

It can also be noted that for the simplified case of f^2 = 
SI3 = SI the #3,4,5 expressions are considerably simplified 
and are given by 



#3 — — — Sli shiKX , 



J4,5 



Sli sin kx ± y SI 2 sin kx + 8S7 2 



for ip = 0. 
(22) 



This means that #3 = #1 for ip = 0; however, we also have 



T 2A 

2A 



for ip = , 
for (p = tt/2 . (23) 



Thus, Li 2 = R\, therefore, the peaks arising from 
sinra = i?i — i.e., #1 — will be completely suppressed for 
ip = 0. Morever, #2 will only appear if 72 is considerably 
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larger compared to all other parameters of the system. 
Numerical study presented in the next subsection sug- 
gests that 72 > 10 A max for the roots #1.2 to start show- 
ing up. It can also be seen that 8\ starts showing up for 
ip = if f^i > 10 f^2 ■ These features can be understood 
by observing Eq. Ijl5(l and comparing the weighting coef- 
ficients of various roots of the numerator and denomina- 
tor. This features are confirmed by the numerical study 
presented in the next subsection. 



B. Numerical considerations 

In the discussion to follow we plots the roots #1,2, 3,4, 5 
as a function of kx and show their connection with the 
behavior of x" vs the probe detuning along the cavity 
field. 

To make contact with our earlier work, ^f, we first 
consider the parameter range with 72 = and = Q3 
and study the effect of increasing 72 on that result. The 
findings are summarized in Fig. In the first column 
we plot the roots #1,2,3,4,5 so that their relation to the 
probe absorption maxima can be established. Then the 
contour-density plot of x" vs the probe-field detuning 
A and the kx and x" vs KX f° r chosen value of A are 
plotted for different values of 72 starting with 72 = 0. 
The color of the line plots corresponds to the horizon- 
tal lines in first-column plots for the respective phase 
value. This correspondence helps to determine the posi- 
tions and number of the peaks in the line plots from the 
places at which the horizontal line intersects the roots 
#1,2,3,4,5(^2). It is clear that the roots $3,4,5, denoted by 
solid lines, are dominant most of the times as opposed to 
<5i 2 ■ The same conclusion can be drawn from the dressed- 
states approach but for different reasons as discussed in 
the next subsection. 

We label the dominant roots, £3,4.5, such that the root 
crossing the A = line (for ip = 0) or the A = line itself 
(for tp = 7r/2) as #3, the root above the A = line as #4, 
and the one below the A = line as #5 as seen in the 
first-column plots in Fig. Among the non-dominant 
roots — denoted by dashed lines in the first-column plots 
in Fig. [3] — the root that coincides with 83 is denoted as 
81 and the other one is 82- It can be noted that the 5x 2 
are independent of the relative phase of the drive fields 
tp. This labeling of the roots will be used for the rest of 
the discussion. 

In Fig. we further observe that with increasing 72, 
the roots £4.5 start diminishing. It is, however, to be 
noted that this behavior can only be seen from the den- 
sity plots of x" an d the plots of #3,4,5 themselves do not 
give this information. Another way to explain the peak 
widths and their dominance is through the decay rates 
of the dressed-states. We discuss the implications in the 
next subsection where we evaluate the dressed-states. 

We consider the results depicted in the line-plots in 
Fig. |3 in further detail. The effect of increasing 72 can 
be easily seen from the peaks arising due to the roots #4,5, 



as seen in the blue plots in Fig. |3] for both the values of 
ip = 0, 7r/2. We first consider the results for ip = — 
For 72 = 0, out of the four blue peaks (corresponding 
to A = 137) occurring the first half-wavelength region, 
the outer ones arise from 84 and the inner ones from S3. 
Thus, the expectation-from the density plots-would be 
that the inner roots would remain sharp and dominant 
while the outer ones will loose their height and sharp- 
ness; this expectation is confirmed by the line-plots for 
72 = 7 and 72 = IO7. The green (A = 57) peaks which 
arise solely through S3 are unaffected by increasing 72. 
The same is true for the red line-plots which correspond 
to the probe detuning of A = giving rise to peaks at 
the nodes of the cavity standing- wave field. Now we con- 
sider the case of ip = tt/2 — Here 83 coincides with the 
zero line, hence for 72 = the red line-plot (A = 0) gives 
equal absorption at all spacial points but starts show- 
ing spacial dependence as 72 increases which can also be 
clearly seen from the density plots in the A = region. 
The blue (A = I67) peaks, in this case, arise from 84 and 
therefore diminish in height as 72 increases. The green 
(A = I27) plots show that #1,2 do not contribute the this 
particular choice of parameters and show zero absorption 
for all values of 72. The results depicted in the line plots 
coincide very well with the corresponding density plots. 

The red curves in Fig. [21 corresponding to the detuning 
A — 0, show different behavior for different phase values. 
For ip — 0, the height and width of the peaks observed 
at the zero detuning of the probe field are immaterial of 
the lifetime of level |2), as they arise from the zero eigen- 
value of the dressed state which does not contain any 
1 2) component. This corresponds to a regime of local- 
ization that is very common in several other localization 
proposals, namely, observance of localization peaks at the 
nodes of the standing-wave cavity field. For ip — tt/2, the 
A = value is special as it does no show any localization 
for 72 = (observe the red plots in the lower row of the 
Fig- 13- It can be seen that these peaks become sharper 
with increasing 72 , whereas the green and blue plots still 
show dimishing height and increasing width of the peaks. 
This can be explained as follows: at A = and 72 7^ 
roots #1,2 dominate as opposed to #3,4,5 for all other val- 
ues of the detunings. Thus, for ip = tt/2, the root 8 = 
starts losing its significance as 72 increases and only the 
nodal points show peaks which arise from #1.2 shown by 
dashed-line plots in first column. This feature is absent 
for (p = as #1 does not occur for the parameters of 
Fig. 0— as explained in the context of Eq. l|2"5)l — and 72 
is not large enough for 82 to show up compared to the 
dominant root #3. However, with a different parameter 
range of values this competition of roots can be seen for 
p = as depicted in Fig. 0] It is also interesting to note 
that for (p — 0, in the current figure, the localization 
peaks at the nodes of the cavity standing-wave field are 
much sharper than the ones observed for p = tt/2. 

Thus, the general conclusion that can be drawn from 
Fig. |31 is that as 72 increases only the maxima due to 
the root #3 (for tp = 0, it) and root #1,2 (for ip = tt/2) 
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72 = 72 = 7 72 = 107 

S X" X" X" X" X" X" 




FIG. 3: The effect of 72 on the locahzation. The parameters are Qi = 30, Q2 = ^3 = 20 71 = 1. Top row: ip — and Bottom 
row: = tt/2. The first coiumn shows the piots of the roots <5i,2 using dashed lines and that of ($3,4,5 by solid lines. The 
colored horizontal lines correspond to values of detuning A chosen to plot x" vs KX m the line-plots shown later in each row 
for different values of 72 shown on top. The line plots are preceded by contour plots of x" to give an idea of its dependence on 
A as well as position, kx. The brightness of a given location in the contour plot is proportional to its height in a 3D plot of x" 
vs A and kx. It can observed that for Q.2 = £1?, = Q the roots 8\ and $3 coincide. The roots £4,5 loose their significance as 72 
increases. This result can be explained through the dressed states approach as described in the text. For A = and 72 = IO7 
the root 8\ — S3 dominates as opposed to 84,5 as expected. Despite expectation the root 82 never appears for ip different from 
7r/2. It can be noticed that for A = and non-zero 72 one might expect both the roots <5i,2 to dominate, however, this is case 
only when ip = n/2. 



show sharp peaks and the other peaks diminish in mag- 
nitude and sharpness. The drive field parameters chosen 
in Fig. |3 were as that of the earlier work in J? except for 
the non-zero 72. 

Now we study the effect of varying the amplitudes of 
the drive fields and go beyond the condition O2 = Q3 on 
the probe field absorption. The results are summarized 
in Fig. 21 The results for = tt/2 in Fig.0]are very sim- 
ilar to the one in Fig. [3 except for the green (A = 127) 
plots. Here the detuning values are the same in both 
the figures, however, the roots £4,5 have a larger range 
than before due to their dependence on the drive-field 
Rabi frequencies. Thus, new roots appear for the green 
plots in Fig. 0] as opposed to no roots in Fig. [3J These 
roots however diminish as 72 is increased making them 
useless for atom localization for larger 72. For the case 
of ip = the roots £3,4,5 have completely different pro- 
files compared to their counterparts in Fig. [3J The red 
plots have the same behavior as in Fig. [3] being a very 
commonly observed localization regime for the probe de- 
tuning A = 0. Moreover, the disappearance of the roots 
£4,5 with increasing 72 exists in this parameter range as 
well. Due to this the blue and green peaks loose their 
height and sharpness with increasing 72. 

We observe that an interesting regime arises when 72 7^ 
and the detuning of the probe field A = 0. In this 
regime the roots S1.2 dominate compared to <53 = for 
ip = tt/2; it can be clearly seen from the last column of 
plots in Fig.m The significance of this is clearly apparent 
for if = tt/2, where the probe absorption is uniform over 
all spatial points for S3 = when 72 = 0; however as £1,2 
become dominant due to increasing 72 absorption peaks 



start emerging at positions corresponding to the nodes of 
the standing-wave field. The behavior for ip = is a bit 
different and it needs complete dressed-states analysis to 
explain that has been found to be quite complicated. 

Another interesting feature observable from Fig. 0] is 
that the root given by S3 is not dominant at all spacial 
positions as it is in Fig. [3J With increasing 72 the rel- 
atively flat regions in the plot of 63 vs kx start loosing 
their significance as 72 is increased. This can be ascribed 
to the broadening of the resonances owing to increased 
72. When S3 remains close close to the line (f2i/2) sinnx, 
the state is very close to the first eigenstate discussed in 
Eq. |^|and these parts remain sharp 72 does not affect the 
sharpness. Whereas, departure of S3 from (Q±/2) sin kx 
lines can be ascribed to increasing components of state 
\a,2), which decays with 72, in the dressed state. Thus, 
the flat regions loose their significance for localization 
with increasing 72. 

Noting that for 72 > A we can expect the behavior 
of the probe absorption to be completely dominated by 
<5i,2 as opposed to £3,4,5 we choose appropriate values for 
the parameters and consider the density plots of x" m 
Fig. |SJ However, this parameter range is unrealistic and 
also not very useful as there will be four peaks observed 
for the detuning lying in the interesting regime as both 
the <5i,2 roots exist giving rise to four intersection points 
in one wavelength for a chosen value of the probe detun- 
ing. Nevertheless, the observations made in the context 
of Eq. Il2i'>[l can be confirmed from the results in Fig. [3] 
Both the curves appearing for ip = tt/2 have equal char- 
acteristics as they both are due to £1,2, which arise from 
the same term in the denominator. However, for ip = 0, 
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72 = 72 = 7 7' = 107 

s X" X" X" X" X" X" 




FIG. 4: Localization characteristics for non-identical drive field intensities and interplay of different roots. The parameters 
are fii = 30, 0,2 ~ 20, O3 = 10 71 = 1. Top row: ip = and bottom row: ip — tt/2. The structure of the figure is the same as 
Fig- CD The roots £4, 5 loose their significance as 72 increases. This result can be explained through the dressed states approach 
as described in the text. For A = and 72 = IO7 the roots <5i,2 dominate as opposed to £3,4,5 as expected for ip — tt/2. This 
can be clearly seen as S3 = line becomes insignificant and only the nodal points (arising from £1,2) remain with increasing 72. 




FIG. 5: Dominance of the roots i?i,2. The parameter val- 
ues are same as in Fig. [3] except for 72. (a) 72 = IO7. (b) 
72 = 10 3 7 (c) 72 = 10 4 7. It can be seen that the less domi- 
nant maxima slowly vanish as 72 is increased. Thus for large 
72 compared to 71 only A = shows probe peak absorp- 
tion. As 72 increases Ri starts dominating which is the same 
as R3, however for somewhat unrealistically larger 72, i?3,4 
roots completely dominate. The apparance of R4 can not be 
explained by the dressed state approach. 



as discussed earlier, Si is cancelled by the roots of the nu- 
merator and instead of Si, the sharper root £3 appears. 
Only when 72 is sufficiently large S3 and S2 acquire same 
sharpness as seen in the last column plots of Fig. |3J 



Deeper understanding of the interplay of different roots 
can be achieved through the dressed states calculation. 
We determine the dressed states in the next subsection 
and explain the above obtained results from a different 
point of view. 

C. Dressed-states approach 



To understand the emergence of several roots for the 
maxima of the probe absorption we consider the dressed 
states approach. The effective Hamiltonian, taking into 
account only the strong drive fields, can be expressed as 



Qgg-iv^ifcscosfla rJi sinra 



Jf eS = - ln 3 e iip e~ ikxcose3 n 2 e [kxcos92 \ (24) 

2 * fii sinra fl 2 e- ikxcos92 



in the basis {\ai) , \a 2 ) , \b}}. Choosing 62 = vr/4 and secular equation 
#3 = 7r/2 + tt/4, in the above Eq. I|24|) we arrive at the 

4A 3 - A(flf sin 2 kx + Vl\ + &1) - flifl 2 ^3 sinracos^ = 0, 

(25) 



9 



where A are the eigenenergies of the Hamiltonian. It can 
be noted that Eq. is identical to Eq. l(2U|l, Thus, 
there is a direct connection between the detuning values 
for the probe field at which it experiences maximum ab- 
sorption, (53,4,5, and the dressed state eigenvalues. The 
dressed state eigenvalues A give the Stark-shifts in the 
energy of the state \a). When this Stark-shifted transi- 
tion |a)-|c) is probed by the weak probe field, the res- 
onances will occur at the points where the probe fre- 
quency matches the energy level difference between the 
Stark-shifted levels \a) and level |c). If the probe field 
frequency, or the detuning A, is chosen such that it is 
in resonance with one of the dressed states then it ex- 
periences absorption maxima. This can be expressed by 
a condition A = <5. It can, however, be noted that only 
the detuning solutions ^3,4,5 can be explained through 
the dressed-states approach and not the solutions (5i,2, 
as explained later. 

Actual form of the dressed states for general parame- 
ters is quite complicated and in not required as we are 
only interested in locating the positions of the resonances 
in the frequency space. In general the probe absorption 
peaks are quite sharp except at the stationary points 
along the x axis and when 72 increases. This calculation 
can be extended further to obtain the spontaneous decay 
rates of the dressed states. These decay rates could then 
give more information about the widths of the probe ab- 
sorption maxima and the loss of sharpness with increas- 
ing 72. 

For a completely general set of parameters evaluat- 
ing the actual form of the dressed states and their spon- 
taneous decay rates is sufficiently involved compared to 
the information that can be gained by such an exercise. 
Therefore, we consider a restricted regime of parameters 
to extract information about the sharpness of the probe 
absorption peaks through the dressed states approach. 
Assuming O2 = O3 = O and ^ = 0we obtain the eigen- 
values to be 

/ — -Oi sin kx , - ^Oi sin kx ± \J8 O 2 + Of sin 2 nx^j | 

(26) 

with the corresponding eigenstates 
where 



with the corresponding eigenstates 



.(±.0) 



e V5 Q 3ftisinfcx ± J8O 2 + f2 2 sin 2 



2Vt 2 + sin 2 kx ± Q,i sin kx J 8 SI 2 + Of sin : 



and jy^' ^ is the appropriate normalization constant. 
Whereas, for ip — tt/2 we obtain the eigenvalues 



W2) 



i|ai) 

— i k x 

— i^e ^2 Oi sin |d2) 



^2 



\b) 

„(±,T/2) 

(±,7T/2) ( ~(±,7T/2) 



0l ) N 

02) 



(30) 



where 



„(±.t/2) 



„(±,t/2) 



fli sin kx ± i\/20 2 + Of sin kx 



i Oi sin kx =F \ 20 2 + Of sin kx 



-20e V2 



(31) 



i Oi sin kx ^ 1/2 O 2 + Of sin 



{0,±-^/2ft 2 + ft 2 sinra} 



(29) 



with JY^ 1 * and tA^''" being the appropriate nor- 
malization constants. 

The message to be taken from the dressed states rep- 
resentation ,jY(c ai \ai) + c Q2 |a 2 ) + Cb\b)) in the bare 
atomic levels is that the decay rate of the correspond- 
ing dressed-state is given by 7 = |c ai | 2 7i + |c a2 | 2 72, as 
the level \b) is the ground state. Therefore, it is clear 
the for ip = the first dressed-state has the deay rate 
7i/4 whereas the other states decay at the rate propor- 
tional to |^ / ( ± '°>| 2 (7i + |4 ±,0) | 2 7 2 ). Resulting in sharp 
localization peaks when the probe field is in resonance 
with the first dressed state and not so sharp localization 
peaks when the probe field is in resonance with the other 
two dressed states. In fact with increasing 72, as seen 
already in the numerical solutions, the latter two states 
contribute wider and wider resonances which are increas- 
ingly useless for atom localization. Similar conclusions 
can be drawn for the case of <p = tt/2; all three roots 
are equally sharp when 72 = and the latter two roots 
increasingly loose their sharpness and decrease in ampli- 
tude for larger 72. This observation can be confirmed 
through the plots in Figs. [21 and 0] 

Another important conclusion that can be drawn from 
the dressed-state eigenvalues is that for the case of ip = 
the eigenvalues can be made to be well separated by 
choosing Qj to be little smaller than Q2 = ^3 = SI. 
In such a case the three roots do not overlap and the 
detuning can be chosen in the range {0, fii/2} to ob- 
tain sub-half-wavelength localization. We illustrate this 
regime in Fig. [5] 

Moreover, this result holds true even when the drive 
fields O2 and O3 do not have the same value, O2 7^ ^3, 
and when fij < 02,0^. Another message to be taken 
from Fig. is that the results for ip — tt are mirror image 
of that of ip — taken around the vertical line kx = 0. 
This holds true for all parameter ranges, hence we have 
plotted only the non-trivial cases ip = and p> = tt/2 
in all the other plots. The range of detuning A spanned 
by the root 63 (for p> — 0,tt) gives an ideal range where 
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FIG. 6: Illustrating the appropriate conditions to obtain good 
sub-half-wavelength probe absorption peaks along the cavity 
field — i.e., sub-half- wavelength localization (observe the green 
plots in last column.). The parameters are fii = 20y, Q.i = 
227, Q3 = 257. 71 = 72 = 7- In the first column we show the 
plots of the roots <5i,2 in dashed lines and that of the roots 
£3,4,5 in solid lines alongwith horizontal lines for chosen values 
of probe detuning A that will be considered for more study 
in column 3. In the second column we show the density plots 
of x" an d in the third column we choose A = 57 (green) 
and show the probe absorption peaks to illustrate the regime 
of sub-half-wavlength localization and its dependence on the 
phase tp — 0. We choose other A values as well shown in 
blue and red to contrast the sub-half-wavelength regime. Note 
f22, 03 > Qi seperates the central root, S3, from the other ones 
#4,5 and provides a phase dependent localization for tp = and 
no localization for tp — n/2. The localization peaks would 
appear in the second sub-half- wavelength region if tp — n. 

sub-half-wavelength localization can be observed, which 
can always be calculated by solving the Eqs. lf%5|) or Q20|l. 
when it is very well separated from the other roots 84^. 
This happens as discussed above when f^i < Q2, ^3. 

In spite of the use of the restricted parameters for the 
evaluation of the dressed-states the results are valid in 
general as our numerical studies show. Neverthelsess, it 



can be noted that the parameter range where 72 has a 
role to play on the dominance of the roots 5\ t 2 or when 
72 ^> S7i,2,3 can not be explained through the dressed- 
states approach (See Fig. [SJ. This breakdown of the 
dressed-states approach for large 72 is easy to under- 
stand. Dressed-state calculation is usually done with 
the assumption that the Drive field Rabi frequencies are 
large compared to all the other parameters of the system, 
which breaks down in the large 72 3> limit, giving rise 
to roots which are not predictable by the dressed-states. 
IV. CONCLUSIONS 

We have studied of a variant of a A-type EIT, where 
a phase dependence is introduced through the applica- 
tion of three driving fields in a loop-configuration. The 
advantage of the phase dependence is in the tunability 
that becomes available to manipulate the response of the 
atomic medium to a weak probe field. By choosing one of 
the drive fields to be a standing-wave field of the cavity 
the phase dependence can be extended to obtain atom lo- 
calization. We have given equations that could be used to 
simulate several, apparently quite different, energy level 
schemes. Effect of different parameters are studied with 
analytical as well as numerical techniques. A dressed- 
states approach is developed and it is used to explain 
the peak probe absorption and the peak widths. Also a 
region of parameters is identified which gives clean sub- 
half- wavelength localization for a wide range of probe de- 
tunings; thus, increasing the applicability of the model. 
In this range of parameters we show how the choice of 
phase governs whether localization would be observed or 
not. 
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